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The amplification of a surface electromagnetic wave by means of ultrarelativistic monoenergetic electron bunch 
running over the flat plasma surface in absence of a magnetic field is studied theoretically. It is shown that 
when the ratio of electron bunch number density to plasma electron number density multiplied by a powered to 
5 relativity factor is much higher than 1, i.e "f rib/n p >> 1, the saturation field of the surface electromagnetic 
wave induced by trapping of bunch electrons gains the magnitude: E x = B y ss 0.16 ulp " lc ( 2 ^ ) 1 / 7 and 
does not approache the surface electromagnetic wave front breakdown threshold in plasma. The surface 
electromagnetic wave saturation energy density in plasma can exceed the electron bunch energy density. 
Here, we discuss the possibility of generation of superpower Teraherz radiation on a basis of such scheme. 
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I. INTRODUCTION 

The surface electromagnetic waves (SEW) on plasma 
surface and plasma-like media (gaseous plasma, dielec- 
tric and conducting media, etc.) attract special atten- 
tion of researchers due to their unique properties. First 
of all, due to its high phase and group velocities close to 
light speed in vacuum at high media conductivity what 
makes them the most valuable in radiophysicsi. The 
SEW are widely applied in physical electronics due to 
its high phase velocity leading to its uncomplicated gen- 
eration by relativistic electron bunches and output from 
plasma. 

Below we discuss the problem of SEW amplification 
with the help of electron bunch running over flat plasma 
surface. We consider the case of ultrarelativistic mo- 
noenergetic electron bunch which remains relativistic in 
the frame of reference of SEW generated by this bunch 
compared to the works^— , where the bunches were non- 
relativistic. Such a problem of generation of three- 
dimensional electromagnetic wave (wakefields) in plasma 
with the help of ultrarelativistic electron and ion bunches 
through Cherenkov resonance radiation was solved u£, 
where it was shown that bunch ultrarclativity influences 



a 'Here we consider the spatially unbounded plasma-(long)bunch 
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significantly the nonlinear stage of plasma-bunch interac- 
tion, in particular, the saturation amplitude of the gen- 
erated wave. 

In the present work we apply the method developed 
ir£ for the case of amplification of a surface electromag- 
netic wave by means of ultrarelativistic monoenergetic 
electron bunch running over the flat plasma surface. The 
interest to the SEW amplification was aroused by its un- 
complicated output from plasma compared to that of the 
three-dimensional wave generated by the bunch as well 
and high magnitudes of SEW energy density. The latter 
is related to the field structure. Thus, as it'll be shown 
below, the SEW saturation energy density can exceed the 
bunch energy density. 

It is noteworthy that the real SEW amplification de- 
vice should be cylindrical what we do comprehend very 
well. However, the problem taking into account the cylin- 
drical geometry is much more complex compared to that 
of plane geometry from the mathematical point of view 
and is not appropriate for illustrative purposes. This is 
why we restrict ourselves to the plane geometry prob- 
lem. Soon, we are planning to finish an article consider- 
ing the real cylindrical SEW bunch-plasma amplifier and 
will present it for publication. 

II. DESCRIPTION OF THE MODEL. DISPERSION 
RELATION 

Let us start our description with the schematic illustra- 
tion of interaction of the ultrarelativistic monoenergetic 
electron bunch with cold isotropic plasma (no thermal 
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motion) being in a rest, which generates the plane wave 
E = Eq exp(—iu>t + ik ■ r), and put the external field as 
absent. @. 

Over the colissionless plasma, filling in the half-plane 
x < 0, with the dielectric permitivity 



Ep = 1 



(1) 



the ultrarelativistic monoenergetic electron bunch, filling 
in the space x > a, with the dielectric permitivity 



£6 = 1- 



<4j 3 

(u> - fc z u) 2 ' 



(2) 



propogates on a distance a. Here u> p = ^4:ne 2 n p /m, 
LUb = yj 4:Tre 2 rib/m are Langmuir plasma electron and 
bunch frequencies respectively (in GSU units) with n p , 
rib being the plasma and bunch number densities in 
the laboratory frame of reference (plasma in a rest) 
(rib *C Up), k z is the longitudinal (directed along the 
velocity of the bunch u) component of the SEW wave 
vector k = (k x ,0,k z ), e= the electron charge, m— its 
mass. The bunch is considered to be an ultrarelativistic 
when 



1 



7 



y/i -U 2 /C 2 



» 1. 



(3) 



The surface wave is a wave of E-type with the nonzero 
field components E x , E z , B y , which satisfy the following 
system of equations^: 



k z 2 E z + %e{x)E z = 
cr 



o 2 e z 

dx 2 

ik z dE z lid . , dE z 



(4) 



where k 2 = k z 2 — elu 2 /c 2 . The system Q is valid for all 
domains shown in Fig. [Tjwith the corresponding substi- 
tutions £ = e p , e = e v = 1, £ = £fe. The electric fields are 
the following functions of the time and the coordinates 



E z = Eq z (x) exp(— iut + ik z z). 



(5) 



Dependence on x is defined by the system ^ and can be 
represented as follows 

Cxe K " x atx<0, 
E 0z (x) = { C 2 e> x + C 3 e~ K » x at < x < a, 

at x > a, 



C 5 e- K » x 



where k p = \Jk 2 - 
n b = \]k 2 - uj 2 e b /c 2 



Sp/c 2 , k v = \/k% — w 2 /c 2 and 



In ultrarelativistic bunches bunch divergence can be neglected 
because when the electron bunch is ejected into the dense plasma, 
within the time t ~ 1/oj p the neutralization of the bunch charge 
occurs prohibitting the bunch divergence 



The boundary conditions can be obtained from the 
field equations by integrating over a thin layer near the 
interface between two corresponding media and have the 
following view: 



E Zp \ x =0 — E Zl/ \ x= o, E Zv \ x=a — E z b\ x =ai 
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Vp\ 



B y J x=0l B 
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In addition to these boundary conditions the following 
condition must be satisfied: 

E Zp \ x — »±oo — E z b\~ >±oo — 0, 



By \x^r±OG B. 



3/£,|->±oo 



0. 



Having solved the system of equations (|4])- ((6]) we can 
finally obtain the following dispersion relation: 



; £ex P (-2oK 1/ )(e6-l)(H--^). (7) 



When the bunch is absent, i.e. rib = and Eb = 1, one 
can get the dispersion relation of surface plasma wave 
from the following equation: 



,y/k* - u 2 /c 2 + ^Jk 2 - u 2 e p /c 2 = 0, 



(8) 



which was studied with the solution ui = loq in detail 
inii. The bunch leads to the amplification of this wave 
and solution of Eq. (JS|) should be found in the following 
form: 



lu = k z ii(l + S) = uj (1 + S), <5<Cl. 



(9) 



Since we took into account that rib <C n p , the highest 
bunch effect on the surface wave occurs when the follow- 
ing Chcrenkov resonance condition is satisfied 



ojq = k z u. 



(10) 



III. ANALYSIS OF THE DISPERSION RELATION. 
BUNCH-ASSISTED GENERATION OF THE SURFACE 
ELECTROMAGNETIC WAVE UNDER THE CHERENKOV 
RESONANCE CONDITION. 

Let us first determine the SEW frequency in a bunch 
absence, i.e. find solution of Eq. ([J| ■ We are interested 
in the frequency range of high-speed waves with luq w 
k z c which can be generated by an ultrarelativistic bunch 
under Chcrenkov resonance condition, i.e. loq w k z u. 
From Eq. (JSJ) follows that such waves can exist only in 
dense plasmas when w « w p and hence e p ~ —lo 2 /ui 2 . 
From Eq. ((8|) we can easily find 
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where the inequality ^ was taken into account. 

Let us now take into account the bunch effect, i.e 
find solution of Eq. ([7]) when the Cherenkov resonance 
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FIG. 1. Schematic illustration of interaction of the ultrarelativistic monoenergetic electron bunch with plasma. Here, the 
corresponding dielectric permitivities are presented as well. 



condition (fTOj) is satisfied. Here, we would like to re- 
strict ourselves to consideration of the ultrarelativistic 
bunch, when \S\~f 2 >> 1, running over and closely to 
the flat plasma surface, 2 3 / 2 awoV / <5/c << 1, for this - 
cxp(— 2a/«„) — I- Then one can get the following solu- 
tion for S 



(12) 



where is the SEW potential and u\ is the bunch elec- 
trons velocity, both measured in the SEW frame of ref- 
erence. In accordance with Lorentz transformations the 
speed of bunch electrons in the chosen frame will be 



ui = -■ 



Ul 



uS'j 



'„ 2 



'„ 2 



L - ^' 7 2 1 - 2(5'( 7 2 - 1) 

-it<5'7 2 « u, at2|<5'| 7 2 << 1 
u/2, at(2|<5'|7 2 » 1), 



The solution of Eq. 
following form 



5 = S'+iS" 



T2l). we are interested in, has the 




It is obvious that the saturation of instability can occur 
when the kinetic energy of electrons, in the SEW frame of 
reference, will become less than the amplitude of the po- 
tential of the plasma wave measured in the same frame. 
In this case the bunch electrons get trapped by the SEW, 
i.e. there will be no relative motion between the bunch 
electrons and the SEW, thus, no energy exchange be- 
tween the bunch and SEW occurs, the bunch and the 
SEW become stationary. For determination of the sat- 
uration amplitude of the potential of the plasma SEW, 
generated by the bunch and amplified with the time with 
the increment 3mj5 = 5" ', we will apply the same method 
used irA Let us choose the SEW frame of reference in 
which the wave is purely potential and its stationary sat- 
uration amplitude $' can be determined by condition of 
the bunch electrons trapping in the wave field^. 
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here the real part of S (5') is considered. 

In the laboratory frame of reference the potentials $o 
and A z arc not zero and 



$o = $'7, A z = -$ - 7$'- 

c 



(15) 



Knowing $o and A z we can determine the fields in the 
laboratory frame of reference: 



E 



_ 1 dA z _ d±a 
c dt dz 
E — -d±o n 



dx 



dA = ^ E 
dx x ' 



>A Z = -ik zl V{l - ^), 
(16) 



It is obvious that the amplified by the bunch SEW wave 
can be easily radiated out of the vacuum domain, in 
which the bunch is running (see Fig. [2]). At the same 
time, the transverse fields E x and B y are much higher 
than the longitudinal fields. This is why we restrict our- 
selves to calculation of the transverse fields in the vacuum 
domain, also because these components form the much 
higher longitudinal radiating component of Poynting vec- 
tor (energy density flux) compared to the transverse one. 



P 



C C 2 

— E r Bi, = — E r 
4tt v 4tt 



(17) 



(14) it is quite easy to calculate the fields E x and B y in the 
plasma domain from Eq. (|16[) . It is noteworthy that 
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in our previous work^ the electric field generated by the 
ultrarelativistic bunch in the range \5\-f 2 » 1 was de- 
creasing with an increase of 7 and we had to determine 
the 7 at which the field is the highest. In this case, the 
behaviour is similar: they increase with 7 at |c>|7 2 << I 
and decrease in the range |<5|7 2 >> I. In the consid- 
ered case the ultrarelativistic range \S\j 2 >> I is of high 
interest. In this range we can obtain the following fields: 



\E x \ = \B y \ = \-^\ = \-iK^ \ 
E^*' = 0.16^(^)V 7 



ui 
u 



(18) 
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These fields are lower than the maximum stationary (sat- 
uration) plasma wave field obtained in£ being equal to 
E m ax = V^ijJpmCy/^f/e. It is assumed that at this mag- 
nitude the surface electromagnetic wave front breakdown 
occurs. In the considered case the transverse SEW is gen- 
erated and its behaviour is not trivial. The only thing we 
can say is that the fields (fT8|) are determined from con- 
dition of the bunch electrons trapping in the wave field 
and that they are lower than those obtained from SEW 
front breakdown threshold. 

Finally, let us determine the SEW saturation energy 
density with the corresponding amplitudes (|I8p 



W 



;(E X 



-B y l ) = —E x 2 = 0.024mc 2 



2n b np /2 \ 



(19) 

These magnitudes can exceed the bunch electrons energy 
density Eb ~ mc 2 ^rib and this fact should not surprise a 
reader (see the Results and Discussions). 



IV. RESULTS AND DISCUSSIONS 

Let us begin with the concluding words of the last sec- 
tion: The saturation energy density of surface electro- 
magnetic wave (|19p . generated by means of ultrarelativis- 
tic electron bunch running over the flat plasma surface, 
can exceed the bunch electrons energy density and this 
fact should not surprise a reader. This is due to the fields 
structure, i. e. |J5 X | = \B y \ » \E Z \. The point is that 
the bunch gets trapped by SEW weak longitudinal field 
component E z , whereas the bunch energy is transfered to 
the whole SEW, i.e. the transverse components E x and 
By get amplified as well forming the longitudinal energy 
density flux. As a result, the SEW saturation energy den- 
sity becomes considerably higher than the corresponding 
bunch energy density. 

The schematic view of the real SEW amplifier is pre- 
sented in Fig. 2. Obviously, the accelerator represents a 
plasma cylinder (plasma is generated in the glass cylin- 
der with the given gas pressure) of fixed length L which 
is blowcd around by the ultrarelativistic bunch. In the 
plasma cylinder of radius r p < (bunch radius) plasma 



of given number density n p is generated. The ultarrela- 
tivstic bunch of radius r& » A/, (bunch width) and of 
given number density rib, the current Jb = 27rr(,A{,ent>it, 
penetrates the metallic vacuumed cylindrical chamber of 
radius R and length L from the left end and comes out 
at the right end being detected by a bunch detector. To 
the right end of the plasma chamber a metallic coax- 
ial chamber is docked where into the SEW, transformed 
into the coaxial transverse electric and magnetic mode, 
is let out. At the junction a partial SEW (longitudinal 
component) reflection occurs and a quasistatic wave with 
increasing along Z-axis amplitude gets formed. Only run- 
ning toward Z-axis wave (forward wave) interacts with 
the bunch whereas the backward wave - does not. 

In conclusion let us make some estimations pursuing 
the goal of employment of the presented above model for 
constructing of supcrpowcrful Tcrahcrz radiator (/g ~ 
0.5 • 10 12 Hz, w = 3 • 10 12 s" 1 ). Since uj = ^ v h th en 
the plasma and bunch parameters can be chosen respec- 
tively. Today, the high-current accelerators are the linear 
accelerators with electron energy of 50 MeV (7 = 100) 
and current density of 500 A/cm" 2 (rib = 10 11 cm -3 , 
net current I = 3 A). The plasma frequency should be 



of order uj v 

n -3 



7^0 



3 • 10 14 s" 1 and 



3 • 10 



19 



cm , atmospheric pressure. Correspondingly, the time 
increment will be ujqS" ~ 2 • 10 8 s" 1 << ojq at a = 0.1 
and the amplification coefficient S"k z ~ 1/L ~ 7.4 • 10~ 3 
enr 1 leading to the system length of 1.34 m where the 
plasma radius is r p = 0.1 cm and beam radius is n, = 0.2 
cm; c/uip — 10 -4 << r p , i.e. the plasma surface can 
be considered as a flat one. Let us notice that the con- 
dition j 2 S' ~ 3 > 1 is satisfied. Finally, let us esti- 
mate the SEW radiation flux. From (TT%|) follows that 
1^1 ~ \B y \ ~ 10 7 V/cm=10 9 V/m, hence, the Poynting 
vector \P\ = c/<ln(E 2 ) ~ 6 • 10 11 W/cm~ 2 . 

Finally, let us notice that presented above results are 
valid only for collisionless plasma, i. c. uq » v e with 
v e being the collisions frequency, satisfying Eq. (fTJ). For 
plasma of electron number density n p and temperature 



T P wc have v e 



10n p 
T. 3 / 2 ' 



where T„ is measured in Kelvin^. 



At cj = 3 • 10 12 s" 1 and 



3 • 10 19 cm" 3 the con- 



dition uio >> v e is satisfied when T e > 2 • 10 5 K. Such 
temperature can be gained during the ionization of a gas 
of atmospheric pressure only with the help of powerful 
lasers. 
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FIG. 2. Schematic drawing of surface electromagnetic wave amplifier. Here, rb » A(, and c/uj p « r p , i. e. the plasma 
surface can be considered as a flat one. 
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